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Abstract 

This paper is devoted to integrable gl(n|m) spin chains which allow for for- 
mulation of the string hypothesis. Considering the thermodynamic limit of such 
spin chains, we derive linear functional equations that symmetrically treat holes 
and particles. The functional equations naturally organize different types of ex- 
citations into a pattern equivalent to the one of Y-system, and, not surprisingly, 
the Y-system can be easily derived from the functional equations. The Y-system 
is known to contain most of the information about the symmetry of the model, 
therefore we map the symmetry knowledge directly to the description of string 
excitations. Our analysis is applicable for highest weight representations which 
for some choice of the Kac-Dynkin diagram have only one nonzero Dynkin label. 
This generalizes known results for the AdS/CFT spectral problem and for the 
Hubbard model. 
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1 Introduction 



Study of integrable systems at finite temperature or at finite volume often lead 
to the functional equations known as the Y-system [1]. This system had been con- 
structed for most of the symmetry algebras with simply laced Dynkin diagrams, see 
for example [2] for a summarizing list. However, a supersymmetric generalization was 
not known until recently when it was proposed for the AdS / CFT integrable system [3] . 
This proposal was then confirmed by explicit thermodynamic Bethe Ansatz (TBA) [4] 
derivation in [5-7]. 

Supersymmetric integrable systems were of interest well before the AdS / CFT case. 
The study of continuous limit of 0l(n|m) rational spin chains first time was systemati- 
cally approached in [8] . Already in this paper a set of TBA equations was derived from 
which one can get in principle the Y-system on a fat-hook (see Fig. 1). Such Y-system 
was however not obtained in [8]. 

In [8] it was understood that, unlikely with a Qi{n) analog, gl(n|m) spin chains 
in covariant representations do not lead in the continuous limit to relativistic field 
theories. One of the reasons for this is absence of the singlet state in tensor product 
of covariant representations, hence the impossibility to introduce antiferromagnetic 
vacuum. Therefore further study in this domain was focused on alternating spin chains 
[. , ]. These spin chains are nonunitary, but as an advantage they do contain the 
singlet state. 

The AdS/CFT integrable system revealed interest to the unitary supersymmetric 
spin chains. Though these spin chains do not posses an antiferromagnetic vacuum, the 
ferromagnetic description is the one that appropriate in the AdS / CFT case. The TBA 
construction leads us to the Y-system defined on a so called T-hook. This is a very 
simply looking system. But first, it should be supplemented with a set of analytical 
conditions [11] which are not that simple, and wc expect to simplify them. Second, 
the derivation of the Y-system is rather lengthy and not at all reflect the apparent 
algebraic structure of the final result. In fact, the Y-system is a kind of the statement 
about the symmetry of the model [12], thus there should be a more direct way for its 
derivation. 

Mainly motivated by the search of the origins of the Y-system, in this paper we 
address the question - how general is the Y-system? For what kind of spin chains do 
we expect to get the structures similar to the one of AdS/CFT case? We state that the 
same kind of Y-systems can be derived for arbitrary gl(n|m) spin chains in a highest 
weight representation and arbitrary choice of the Kac-Dynkin diagram, with condition 
that there is a choice of grading when only one Dynkin label is not zero^ . 

The usual way to derive the Y-system is to go through the whole TBA procedure. 
However we will show that a T-hook pattern can be spotted already at the first step 
of TBA, when one formulates the string hypothesis and writes down the set of linear 
integral equations to approximate the Bethe Ansatz. 

String hypothesis has a long story of application in the Bethe Ansatz systems. It 
already appeared in a seminal work of H. Bethe [13] where it was successfully applied 
to support hypothesis about completeness of the Bethe Ansatz. Later on, it was used 
to study spectrum of excitations of various spin chains around antiferromagnetic vac- 

^If this label is equal to one, such representation is sometimes called fundamental. This 
terminology will not be used in this article to avoid confusion. 
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uum [14,15]. The string hypothesis apphed for AdS/CFT case in [16] originated from 

the work of Takahashi [17] on the Hubbard model, where stack configurations (trian- 
gles in this paper) were first time introduced. The string hypothesis was used in [8] 
when continuous limit of 0l(n|m) spin chain was considered. 

The main claim of this paper is that the integral equations which are obtained in 
the thermodynamic limit and under assumption of string hypothesis can be rewritten 
as in a remarkably symmetric form: 

X] * Po-',s + X] * P*o.,s' = source term. (1) 

a' s' 

These equations cover all the cases of integrable chains with a symmetry group from the 
A series. Up to our knowledge, similar symmetric formulation, in a Fourier space, first 
time appeared in [18] for gl(n) case, here we make its supersymmetric generalization. 
The operator „/ is decomposed to 

Ka,a' = Sa,a'K — Sa,a'+1 " Sa,a'-1, (2) 

where definition of K requires some explanation and is discussed shortly below. 

Pa^s and p* g denote the density functions for particles and their holes^ , and we 
label different types of particles by two integer indices {a, s}. Exact rule for assigning 
{a, s} labels to different particle types is given in the main text. This rule leads to 
organization of excitation types into one of the patterns of Fig. 1, depending on a real 
form of a symmetry algebra. 

Equation (1) is satisfied everywhere inside of one of the domains of Fig. 1. There 
is an additional equation in the corners of fat- and T-hooks which is roughly obtained 
by inversion of Ka^a' and Kg^s' and is given by (61). 

Equation (1) already appeared in the PhD thesis of the author [19], however with 
an indigestible proof and without complete understanding of the structure of equations 
in the corner. Also the relation between different choices of Kac-Dynkin diagrams was 
not clarified. In this paper with the help of the fermionic duality transformations we 
resolve all the complications that appeared in [19]. 

To define function {K * p) {u) we first introduce the resolvent R of p: 

oo 

R{u)= [ dv^^, -J-{R{u + iO) - R{u-iO)) = p{u), (3) 

J U — V ZTTZ 

— oo 

and functions 

}C±R={D + D-^)R, Im(M) ^ 1/2, (4) 

where D = e*/^^" is a shift operator. The functions /C+i? and IC-R are defined in the 
domains Im(u) > 1/2 and Im(u) < — 1/2 respectively and analytical there; then they 
can be continued to the whole complex plane, the result may depend on the chosen 
path. K * p is defined by a difference of K.+R and K—R: 



^In the horizontal strips of Fig. 1 p is assigned to particles, and p* - to holes, while the 
assignment in the vertical strip is reversed. In the central region assignment rule depends on 
a particular choice of a Kac-Dynkin diagram. 
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a) 5u(n), strip 



b) 5u(n, m) , slim hook 




c) su(n|m), fat hook d) su(n|m|p), T-hook 



Figure 1: Positions of different types of particles (string configurations) that 
define equation (1). The same domains appear for the Y-system (8) and for the 
Hirota relations (9) between the transfer matrices. 



K*p = — ^—{IC+R-K-R). (5) 
2iTi 

If p{u) is an analytic function in the strip |Im(u)| < 1/2 than we get 

{K*p){u)=p(u+'^+p{u-'^ (6) 

for real values of u. However, if p{u) has branch points in this strip than analytical 
continuation of 1C±R is not unique and K * p would be a multivalued function. 

As was discussed in [19], from equations (1) the y-system can be straightforwardly 
written: 

1 + =(l + l^a,.')^-'*, (7) 

where Ya^s = is evaluated for the configuration that minimizes the free energy. 

If Ya^s is analytic in the strip |Im(ti)| < 1/2, which is so for relativistic models,^ 
than (7) reduces to a standard formulation of the Y system^'^: 

Y+Y- - (1 + Wl)(l + ^a,s-l) .g^ 



^Though in this paper we consider spin chains, relativistic models for which y-system is 
written may be obtained as a continuous limit of spin chains. Thus the equations (1) are the 
same for relativistic models after we exchange p and p* for momentum-carrying excitations. 

^notation f^{u) — f{u ± i/2) is used. 

^We were informed that independently V. Kazakov, A. Kozak, and P. Vieira obtained Y- 
system (8) for a general T-hook pattern. Due to relation between (7) and (1) this is the same 
observation we discuss in this article. For su(2, 2|4) case this is of course the observation of [3]. 
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For the AdS/CFT case some of F- functions have cuts in the mentioned strip or on 
its border, therefore equations (8) are satisfied only for the particular choice^ of ana- 
lytical continuation of functions IC±R. A possibility to choose different prescriptions 
for continuation is a source of seeming discrepancy between works [7,20] and [3,6]. 
Let us note that all possible pathes should be chosen in order to define completely the 
Y-systcm while (8) corresponds only to one of the pathes. Other choices should lead 
to monodromy constraints on the y-functions formulated in [11]. 

Domains of Fig. 1, obtained here through the analytical study of sting hypothesis, 

have a quite different algebraic interpretation. One can assign so called rectangular 
representations (see sections 2.4 and 4) to each node of these domains in a way that 
transfer matrices in these representations would satisfy Hirota bilinear equations: 

Ta,sTa,s = Ta,s+lTa,s-l + Ta+l,sTa-l,s- (9) 

The construction is well defined in the case of spin chains. In the case of relativistic 
models it was built explicitly only for some particular examples [21], and this allowed 
to find spectrum of the excited states''. 

Besides rectangular representations, there is a way to define the generalization of 
the Young tableaux that can be inscribed into a T-hook such that it defines a unitary 
highest weight representation (UHW) of the symmetry algebra [12]. In this article we 
argue that all UHW can be obtained in this way. 

The article is organized as follows. In section 2 we introduce notations for gl(n|m) 

algebra, define the Bethe Ansatz equations and recall realization of fermionic duality 
transformations on the Bethe Ansatz equations. Section 3 contains the derivation of 
our main result - equation (1). In section 4 we explain how to parameterize UHW 
in terms of generalized Young tableaux and conjecture that this allows to describe 
all UHW. We discuss obtained results and outline possible directions for future study 
in section 5. Article contains also three appendices, first two are devoted to showing 
some evidence supporting the string hypothesis. The third one shows how to treat 
more involved situations with the source term. 

2 Integrable supersymmetric spin chains 

This section collects known in the literature results about Bethe Ansatz for rational 
spin chains. Since in this paper we are dealing with the string hypothesis, we will 
need basic knowledge about representations which allow for this hypothesis. These are 
so called rectangular representations, name coming from the shape of partition that 
encodes them. 

^for some small a and s additional kernels that depend on the 't Hooft coupling constant 
should be added to Kay and Ks^s> . However, for the particular choice of analytical continuation 
this does not affect (8). 

^Ideology of TBA [4] allows for computation of ground state only, though there are analyt- 
ical continuation tricks [22] which allow computing energies of some excited states. Another 
approach is to derive nonlinear integral equations starting from the lattice regularization [23,24] 
which can be generalized to compute energies of excited states [25]. 



6 



2.1 Conventions for 0((n|m) algebra. 

The 0[(n|m) algebra is described using its generators Eij, 1 < i,j < n + m, which obey 
the supercommutation relations: 

[Eij,Eki] = 8,kEa - (10) 

It is useful to think about E^j as the basis vector in the space V ®V* , where V = C""*""^ 
is a graded vector space. Than \i\ means the grading of the basis vectors Vi in V and 
can be (even) or 1 (odd). 

There are different ways to distribute gradings in V . We are not going to make 
some distinguished choice. Instead, we use the Kac-Dynkin diagram to define the 
grading. If the z-th node of the Kac-Dynkin diagram is white then vi and Uj+i have 
the same grading, if this node is crossed - then Uj and Ui+i have the opposite grading. 
Together with agreement what is the grading of vi, this uniquely defines Then, the 
generators Eij are odd if \i\ ^ \j\ and even otherwise. 

Once the basis vectors Vi are chosen, we define Borel decomposition of the 0l(n|m) 
algebra be requiring that E'^j with i < j are raising operators (correspond to pos- 
itive roots). This will be called a standard Borel decomposition with respect to a 
corresponding grading. 

The irreducible highest weight representation is uniquely defined by 

Eij\n) = 0, i<j, 

Eii\n) = mi\n). (11) 

To denote the weight of the irrep we will use the square brackets: [mi, . . . ,miv+M]- 
To render formulas simple, the hatted notations will be also used: 

Eii = {-l)\'\Eii, fhi = {-l)\'\mi. (12) 

The general rule is that we will put hats when factors of type (—1)'*' are used in the 
definition and no hats otherwise. 

The Cartan generators hi of the s[(n, m) subalgebra and correspondingly the Dynkin 
labels uji are defined by^ 

hi = (— l)'*'^j = -E'm — , ^^^\ 
Qi = (-l)l*la;i = rhi- fhi+i . 

Set of the Dynkin labels will be enclosed in the angle brackets: (cDi, . . . jcZj^v+M-i)- 

Note that negative value of cDj would still mean compactness if|i| = |i + l| = 1. 
And inverse, if \i\ = |i + 1| = 1 and cDj > 0, the corresponding representation is infinite 
dimensional. 

2.2 Bethe Ansatz equations 

Solution of the integrable spin chain by the inverse scattering method is based on the 
procedure of diagonalization of the transfer matrix. The latter can be diagonalized by 
means of the (nested) algebraic Bethe Ansatz, once each site of the spin chain is in the 
highest weight representation. 

^Sorne of the Dynkin labels may differ by sign from those defined in [26]. 
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The eigenvectors are parameterized by the solutions of the Bethe Ansatz equations. 
These equations for a generic g[(n|m) spin chain are given by [27,28]: 

n+m-l Mfc ^(t) (k) 



li li {<?) (fe) i, V . / W i ^ • 

k=l j = l Ui -Uj ~ 2\^^^ "fc ) ^^+lK- - 



Here £ = l,n + m— 1, i = l, Mi and the following notation is introduced 
• Ai is defined by 

L 
k=l 



Mu) = e'^^l[iu-Xk + ifhl^^)- (15) 



• L is the length of spin chain. 

• A grading of the gl(n|m) and its Borel subalgebra are chosen, are the simple 
positive roots with respect to this choice. The scalar product between simple 
roots (graded Cartan matrix) is given by: 

( a,, ak ) = Se,k{{-lf^ + (-1)1^+^1) - (-l)l^+^l%i,fe - (-l)l^l5,_i,fc. (16) 

• Mk is the number of Bethe roots at the Dynkin node k. 

• [rrij^\- ■ ■ , 'fn^k'^^^] is the weight of the highest weight vector at spin chain cite k. 

• C( are defined by 

c,=x:(-i)i'i- (17) 

k=l 

• Xk c^i'G the inhomogeneity parameters - arbitrary numbers which are assigned to 
each node of the spin chain. 

• e*'^^ are the twist parameters. They appear if we insert in the supertrace which 
defines the transfer matrix the term Q = diag(e*'^i, . . . , 6*'^'^+'"), where e^'^^ = 

Twists and inhomogeneities are auxiliary parameters. Variation of them allows 
keeping solution of the Bethe Ansatz equations in a general position. We will 
use dependence on twist for study of string hypothesis in appendices A and B. 

The weight [fii,. . . , ^in+m] of the eigenvector associated with the solution of the 
Bethe Ansatz equations is given by 



lii = Mi_i-Mi+Y,m^^- (18) 
k=l 



The eigenvector is the highest weight state in the case when (/)£ = 0. Adding Bethe 
roots at infinity corresponds to the action by the lowering operator. When the twists 
are turned on, the Bethe roots at infinity become finite^. 



^When the twists are present, notion of highest weight is not applicable since the symmetry 
algebra is reduced to the Cartan subalgebra only. 
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rh' 



m 



a+b+1 



Figure 2: Fermionic duality transformation for the case 6 7^ 0. In the case b — 
a, b, c do not change. Horizontal direction corresponds to generators with 
|z| = 0, vertical - with |z| = 1. 



2.3 Duality transformations 

Consider change of the basis in V: 



(19) 



where (Tk{k) = + 1, (Tk{k + 1) = A; and C7fc(i) = i for i ^ k,k — 1. 

This transformation changes the standard Borel subalgebra and therefore changes 
notion of the highest weight representation. So, operator e+ = Ek^^^i transform to 
e_ = Ek+i^k- The highest weight vector |0) remains such only in the case if 



hk\n) = o. 



(20) 



In all other cases \Q) does not remain the highest weight, however representation itself 
may remain of the highest weight type. Let us consider different possibilities 



Compact case. If \k\ = |A; + 1| and is a nonnegative integer then su(2) subal- 
gebra, generated by e± and h^, forms a finite dimensional representation by acting on 
It is straightforward to see that the vector e^''\Q,) becomes a highest weight vector 
with respect to the standard Borel subalgebra of (Tfe(gt(n|m)). The weight of e^''\Q) in 
correspondence with ak{h) is the same as the weight of in correspondence with h 
(here h is the Cartan subalgebra). 



Fermionic case. If e± are fermionic generators {\k\ / |A; + 1|) then together with 
hk they form an su(l|l) subalgebra. By acting on it forms two-dimensional repre- 
sentation. And it is straightforward to check that e_|r2) is the highest weight vector 
with respect to the standard Borel decomposition of (Tfc(g[(n|m)). Transformation of 
the highest weight, assuming that = and + 1| = 1, is given by 

[. . . , fhk-ufhk, fhk+ufhk+2, . . .] ^ [. . . , fhk-urhk+i -%mk- 1, mk+2, • • •]• (21) 
On the level of the Kac-Dynkin diagram the transformation is shown in Fig. 2 



Noncompact case. If = j/c + 1| and ujk < or non-integer, the highest weight 
representation of the su{2) subalgebra is infinite dimensional and becomes of lowest 
weight type after ak transformation. For the whole 0[(n|m) algebra, representation. 
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which is of highest weight type for the standard Borel decomposition with respect to 
initial grading, is no more of highest or lowest weight type after transformation 
applied to a noncompact bosonic node. 



Let us now consider how the automorphism is realized on the level of the Bethe 
Ansatz equations. 

Fermionic duality. First consider the case when e± are fermionic. The transfor- 
mation known as a fermionic duality transformation was studied number of times in 
the literature [20-32]. 

Since we heavily rely on the fermionic duality transformation in this text, we will 
repeat basic facts about it. Assume the grading as is in the l.h.s. of Fig. 2, labeling of 
Bethe roots as in (48) and consider Bethe equations for u, 9, v nodes. To be definite we 
take u, V to be bosonic nodes, though this is not of importance. The Bethe equations 
are written as 

n- (0++0r A + 

(22a) 
(22b) 
(22c) 







K^2 


Qu- 








Qu Q~v 


An 




Qu Qv 


Ae' 


Q\ 


IQv Qv2 


_ 


Q\ 


rj Qv^Qv2 


Ai ■ 



Here Qa are the Baxter polynomials, the first equation is evaluated at zeroes of Qu, 
the second - oi Qo, and the third - of Q^,. 

It is possible to simultaneously translate all the Bethe vectors by a constant. This 
allows changing all q by a constant. We use this freedom to put cg = 0. Therefore 
Cu = Cy = This explains the shifts by i/2 in A in first and third equations of (22). 

Equation (22b) can be rewritten QQ relation: 

AeQ-Qt + KQtQv = QeQe (e^^^ + e*'^") • (23) 

The requirement that Q-g is a polynomial is equivalent to the Bethe Ansatz equations 
(22b). Note now that Qq and Qq enter in a symmetric way in (23), thus a natural idea 
is to exchange the role of them. 

First, one can write the Bethe equations for the zeroes of Qq. These equations are 
exactly the same as (22b). To preserve the pattern of (14) we write it as 

91% = (24b) 
QuQ 

This equation is evaluated at zeroes of Qq - dual fermionic roots. 

Let us replace, using (23), Qg with everywhere in (22a) and (22c) . For each 
application of (23) when replacing, only one of the terms in the l.h.s. of (23) survives 
and the result is the following: 

Qu-2Qq A/, 
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Equations (24) precisely match (14) for the choice of weights given by the r.h.s. of 
(21). 

If Kq = A„, the A-s cancel out from (22b) and therefore r.h.s. of (22) remains 
unaffected after duality transformations. This corresponds to the case hk\^) = in 
which \Q) remains highest weight vector after duality transformation. 

Sometimes fermionic duality transformation is sought as a particle-hole transfor- 
mation for fermionic Bethe roots. This is not completely true as only part of dual 
fermionic roots are holes, see section 3.3. 

Bosonic duality. It is possible to perform duality transformation in the bosonic 
case also [33-35]. It is based on the fact that Baxter equation allows for two solutions 
which are related by the Wronskian relation. Duality transformation exchanges these 
two solutions. In the compact case both solutions are polynomials and the pattern 
of the Bethe Ansatz equations after duality transformation does not change. In the 
noncompact case second solution is not a polynomial. 

2.4 Rectangular representations 

For general choice of the highest weight representations the Bethe equations are not 
invariant under complex conjugation and therefore do not admit real solutions. The 
equations are however invariant if the weights satisfy one of the following two condi- 
tions^*^: 

either m^^^ = + q (25a) 

or m^^) = m^^+^l (25b) 

Since we are aiming to formulate the string hypothesis, the reality of solutions 
is a necessary requirement, therefore we stick to the representations obeying (25). 
Moreover, we will consider a simpler case of representations - those that admit grading 
in which (25a) is true only for one instance of i. Such representations are characterized 
by the property that only one of the Dynkin labels is nonzero. On the language 
of generalized Young tableaux inscribed into T-hook (see section 4) these are the 
representations given by a rectangle defined by position of its two corners: one at 
{0,0}, another at {a,s}^^. Due to this such representations are called rectangular. 

Rectangular representations are all unitarizable, and as a consequence of unita- 
rizability, their tensor power decomposes into irreps. For gl(n|m) case this is not a 

^°To be strict, invariance under complex conjugation may be achieved also if there are 
special relations between weights of different nodes of spin chain. There is an example in the 
literature [9] which fits into this class. Interestingly, solutions of Bethe Ansatz equations in [9] 
turned out to be not always real. We consider only spin chains with the same representation 
at each node. 

"'^"'^If {a, s} corner of the rectangle belongs to right or left band of the T-hook, it defines the 
representation with more than one nonzero Dynkin labels. However, if we imbed the symmetry 
algebra in a larger one such that {a, s} is in the vertical strip, it would be possible to choose 
grading where all except one Dynkin labels are zero. See appendix C for more details. 
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trivial property. In fact, appearance of indecomposables in tensor products of irreps 
in 0[(n|m) is one of the reason that the resuhs for 0[(n) case cannot be directly gener- 
alized. But rectangular representations allow for such generalization, and this what is 
basically done in this paper. 

Another appealing property of rectangular representations is that they allow for a 
simple formulation of the unitary scattering matrix. Indeed, one can choose a spec- 
tral parameter in such a way that Lax operator satisfies L{6)L{—9) oc I. The latter 
relation is a consequence^^ of the property oc J^, where j\ is the Lie algebra gen- 
erator in a particular representation. The property oc holds for rectangular 
representations, but is not true in generic case. 



3 String hypothesis and functional equations. 

3.1 Status of string hypothesis. 

Let us explain a standard argument for appearance of string configuration on the 
example of the XXX magnet. The Bethe equations are written as 

where L is a length of the spin chain and M is a number of Bethe roots. 

The argument is the following. Suppose that L is large and some Bethe root n„ 
has a positive imaginary part: Im(n„) > 0. Then the l.h.s of (26) is exponentially 
large with L. To achieve this large value on the r.h.s. there should be another Bethe 
root Un' ^ Un — i, with the help of which the pole in the r.h.s. is created. Repeating 
the same arguments for and using the reality of solution of the Bethe Ansatz [36] , 
we conclude that the Bethe roots are organized in the complexes of the type: 

Uk = uo + ik, k = 2~''"'~2~' ^'^'^^ 

where s is an integer. These complexes are called s-strings. 

String hypothesis in its strong form states that all solutions of the Bethe Ansatz 
equations can be represented as a collection of strings, and that are approximated 
by uq + ik values with exponential in L precision. 

In its strong form the string hypothesis is wrong. However there is an evidence that 
its weaker version is correct if the proper thermodynamic limit is taken. The weaker 
version states that most of the Bethe roots are organized into strings with exponential 
in L precision, and that the fraction of solutions which significantly differ from (27) 
decreases to when L — )• oo. We discuss in more details applicability of the string 
hypothesis in appendix A. 

If the string hypothesis is valid at least in its weak form, and if energy of most 
of string-like excitations is low^^, the approximation by strings is a good one. In 
the following we will work using this approximation. This allows us to rewrite the 

^^We thank to Carlo Meneghelli who pointed out this explanation to us. 
^^this is so for the antiferromagnetic Hamiltonian. 



12 



Bethe Ansatz equations into equations for centers of strings (called string Bethe Ansatz 
equations henceforth). 

Further study, as it will be clear below, shows that approximation by string Bethe 
Ansatz equations is good when Bethe Ansatz can be approximated by the linear integral 
equations for the densities of Bethe roots. We are going to derive now the functional 
equations equivalent to the linear integral equations, and the final result will be given 
by (!)• 

3.2 Holomorphic projection. 
3.2.1 su(2) case. 

In the large volume limit the Bethe Ansatz equations can be approximated by a func- 
tional equation. Let us show first how this is done for the solution of the 5u(2) spin 
chain (we are following [19] where more details can be found). The starting point is 
the Baxter equation 

W+Q— + W-Q++ = TQ, (28) 

M 

where W = u^, Q = Yl{u — u^), and T is the transfer matrix. 
k=l 

The ratio r = ^'^'^^^ is assumed to have the property for L ^ 1: 

|r| » 1, Im(it) > 0, 

|r| < 1, lm(u)<0. (29) 

This assumption is obviously not valid for large values of u. Note however that string 
hypothesis is also not valid at large values of u. For u closer to origin (n < \/L) 
assumption (29) is a good approximation. Moreover, below we will show how to derive 
exact string configurations from (29). 

Based on (29), we can simplify Baxter equation to 

W+Q— = T Q, Im(M)>0, (30a) 
W-Q++ = T Q, lm(«)<0. (30b) 

We call this approximation the holomorphic projection of the Baxter equation. We see 
that the l.h.s. of (30a) has zeroes at position u^ + i where is a Bethe root. R.h.s of 
(30a) should posses the same zeroes, therefore each such zero should be either zero of 
Q or zero of T. 

Accompanying roots. Suppose first that all Bethe roots are real. Therefore u^ + i 
cannot be zero of Q, thus it is a zero of T. Making the same consideration for equation 
(30b) we conclude that T has also zeroes at u^ — i. We call u^rti accompanying zeroes. 
Therefore T can be represented as 

T = Q++Q--Q*. (31) 

The zeroes of the polynomial Q* = W^i{u — Uh^i) are called holes. All the holes should 
be real, otherwise Uh^i ± i will be the Bethe roots in virtue of (30a) and (30b). 
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By introducing the resolvents 

R = dulogQ, R* = du\ogQ* (32) 

and taking the logarithmic derivative of (30a) and (30b) we get 

(1 + D+2)i? + i?* =D+'^^, lm(u)>0, (33a) 
{1 + D-'^)R + R* =D-^^, lm(u)<0. (33b) 



These equations are equivalent to the linear integral equation which describe XXX 
spin chain in the thermodynamic limit (see introduction to part III of [19]). 

Including string configurations. If we allow the Bethe roots to be complex 
then the roots it i may be the Bethe roots as well. In theirs turn the Bethe roots 
Uk^i may imply presence of ± 2i roots e.t.c. In this way we see that Bethe roots 
are organized in string configurations once (29) is satisfied. 
Let us rewrite the Baxter polynomial in the following form: 

oo Ms 
s=l k=l 

where Uk,s is the center of s-string. Here we introduced the "D-number" 

[s]d = = D-'+^ + D-'+^ + . . . + D'-^ (35) 

and used notation 

QHD)=^HD)logQ^ (36) 

Repeating the logic of the accompanying roots, we conclude that the holomorphic 
projection of the Baxter equation leads to 

oo 

W+ = Q*l[Qf'^'^^'-\ (37) 

s=l 

The latter equation is written in the upper half plane. Equation in the lower half plane 
is just a complex conjugated one. We keep agreement to write only u.h.p. equations 
in the following. 

The functional equation, the logarithmic derivative of (37), is given by 

oo J- 

{D + D-^)J2d'Rs + R*i=D- (38) 

s=l 

with definition of resolvents analogous to (32). 
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Functional equations for string configurations. We obtained equation (38) 

which can be considered as equation for the density^^ of 1-strings in the presence of 
higher-length strings. We will now derive equations for the densities of s-strings. 

For this we should write Bethe Ansatz equations for the center of strings. They are 
obtained by multiplication of the Bethe Ansatz equations for the Bethe roots - centers 
of strings. The result is the following 



J's{uk,s) = 1, (39) 

where we defined 

oo 
s'=l 



-lYT.(u) = W"--"- n qT""'"'"' , (40) 



D-D-^ ^ ' 

Then we introduce a Baxter-like equation, perform its holomorphic projection and use 
the argument for accompanying roots to get functional equations. The net result for 
the rules to obtain functional equations is the following: 

1. Write Bethe Ansatz equation in the form 

7{u,) = \, T{u) = n T,(D-^) . (42) 

k fk 

where fk is some function of u (in the case of (39) this is or Q^), and 7^(t) is 
a polynomial (or more generally, series) in t such that 7^(0) = 0. 

2. The functional equations for the resolvents Rk = 5„log/fe and R = Y1 
then written in the form: 

R + R* = ±^V{D+^)Rk, lm(u)>0, (43) 

k 

and the conjugated equation in the l.h.p. 

3. The sign is chosen to be positive if |J^(tt)| S> 1 for Im{u) > and to be negative 

if \J='{u)\ <C 1 for Im(n) > 0. 

Application of these rules to (39) gives: 

oo J- 

Rt + {D + D'') ClMs' = D^- . (44) 

s'=l 

Now note that (-C"*")"! = Cgs', where 

Css' = {D + D~^)5ss' - 5s,s'+i - Ss,s'-\, s, s' > 1. (45) 
Therefore we get 

oo J- 

iD + D-')R, + YC,s'R*s = S,,i-. (46) 

s'=l 



*Once the resolvent is known, the density is restored using (3). 
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3.2.2 su(n) ceise. 



Exactly the same procedure can be applied for the Bethe equations in the su(n) case, 
and in general any compact rectangular representation {A, S} and the case when in- 
homogeneities Xk present can be considered. The result is 

n— 1 oo 

Caa'Ra',s + ^1 ^ss'R*a,s' = ^a,ASs,sJ , (47) 
a'=l s'=l 

where Ra^g is the resolvent for the density of s-strings composed from the Bethe roots 
that corresponds to the simple root of the Dynkin diagram. J stands for the 

L 

resolvent of a source term and is given by J = H u-v ■ 

k=i 

If we take the difference of (47), defined for Im(u) > 1/2, and the complex conju- 
gation of (47), defined for lm{u) < 1/2, we get exactly the equation (1). To define this 
difference, a prescription for analytical continuation is needed, but this is also the case 
for (1). 



3.3 String hypothesis for g((n|m) case. 

In the supersymmetric case a new type of the Bethe equations is added - equation for 
the fermionic node. Assume the following labeling of the Bethe roots: 



-<y<y<ymy<yo- 

■ ■ ■ U-2U-1 Q Vi V2 • ■ ■ ' 

and for simplicity n_i = u, vi = v. 

The number of Bethe roots decreases from left to right: . . . > > Mg > > 
. . .; this ordering defines what term in the Baxter equations should be suppressed for 
the holomorphic projection. 

The Bethe equations for the 0-roots are written as 

The QQ-relation is the the following: 

QlQv - QZQt = QeQe- (50) 

Here Qu,Qv,Qe stand for the Baxter polynomials of the corresponding Bethe roots 
and Qq is for dual fermionic roots, see section 2.3. 
The holomorphic projection gives 



QtQv = +QeQe, lm(u)>0, 

QZQt = -QeQe, lm(u)<0. ^ ^ 



We see that since there is no term on the l.h.s. of the QQ-relation, fermionic 
Bethe roots cannot form string configurations by themselves. Suppose however that 
among the u-roots there is a string of length s with center at uq. Holomorphic pro- 
jection (51) implies then that there is a string of the length s — 1 (with the same uq 
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center) which should be the string^^ of Q0 or Qg. Equivalently, presence of s-string 
among ti-roots imphes presence of {s + l)-string either in Qq or Qg. 



3.3.1 Triangles (distinguished diagram). 

Let us focus on the ease of s-string of f-roots. Suppose that induced (s + l)-string is 
that of Qg. Then, since 6 roots do not interact with v roots, there is no back reaction 
of (s + l)-string, and s-string of u-roots is just a string configuration discussed in the 
bosonic case. Situation is however different if (s + l)-string is the one of Qg. Let us 
perform the duality transformation on the 9-node. Then i;-node becomes fermionic, 
and it still contains the s-string. Fermionic node cannot contain s-string by itself, and 
there is also no (s-l-l)-string of 9 that would induce the s-string of vs. Therefore there 
should exist an (s — l)-string of V2. The procedure of dualisation is applied again, now 
to the node vi, and so forth until we reach the 1-string at the node Ws-i. Coming back 
to the variables before dualisations, we conclude that an (s -|- l)-string of fermionic 
Bethe roots may exist if it is a part of the following "triangle" configuration: 



(52) 



^8KD a 

Vi Vs-1 



We will call this configuration (s -|- l)-triangle. 



3.3.2 Trapezia (non-distinguished diagram). 

The trick with fermionic duality transformation can be applied to reveal another pos- 
sible configuration. Let us come back to the case when there is an s-string of v-s and 
an induced (s -|- l)-string is that of 9. By performing the duality transformation on 
the 9 node we get the so called trapezium configuration which consists of s-string of 
v-s and (s -|- 1) string of 9. Note that v becomes a fermionic node after the duality 
transformation. 

The general configuration called (s, s')-trapezium is described as follows: 





• 
• 








• 
• 




I. 



(53) 



" ^s'-Fl 

""^^Thc reader may ask why the string cannot be distributed among both Qg and Qg. The 
reason is in a back-reaction of 9 Bethe roots in the Bethe equations for u. This is most clearly 
seen through the duality transformations. 
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This configuration can be obtained after for example by the set of duahty transforma- 
tions dl or d2: 






(54) 




The remarks follow: 

- There are some complications if the duality transformation is applied to the node 
with the source. In general situation, fermionic nodes - boundaries of trapezium 
cannot be separated by the node with the source. 

- Any trapezium becomes a bosonic string for a certain choice of Kac-Dynkin 
diagram. To see this, we just make duality transformations dl or d2 back. We 
have a choice to obtain either s or s' string. 

- Triangle can be always seen as unfinished trapezia. For this we need to enlarge 
the Dynkin diagram (52) with one auxiliary (virtual) fermionic node. The en- 
largement can be always done, with the additional requirement that there is no 
Bethe roots at the virtual node. After this trick, we can make duality trans- 
formations similar to the inverse of d2 to bring the s-triangle to the s-string. 
However, virtual fermionic node will be dualized, and we should require maxi- 
mal excitation at the dualized virtual node. 

- Duality tricks^^ allow us to make even a more general transformation. An s- 
string at node k for one choice of the Kac-Dynkin diagram can be represented 
as an s'-string at node k' for some other choice of the diagram ii s — k = s' — k'. 
To perform this transform we may probably introduce a virtual fermionic node. 

3.3.3 Diamonds (noncompact, both gl{n) and gl(n|m)). 

Consider now noncompact case and the diagram of the form: 



where n is a momentum-carrying node. If the nonzero Dynkin label m is positive, the 
Bethe roots of type u may form s-strings. After dualisation of 6i and 02 nodes the u 
node transforms to noncompact one and the s-strings become a set of two trapezia. 
This shape will be called diamond. The general (s, s', s")-diamond configuration looks 

^^Somehow similar tricks were used in [37] to find a new class of explicit solutions of 1-loop 
Bethe Ansatz for AdS/CFT spectral problem. 




(55) 



$2 u 6*1 • • • ' 
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as follows: 




"i^^ ^^J.'^ , (56) 

^K>©<X>8K> 

(s, s', l)-diamond can terminate also on the bosonic nodes, as triangles do. In 
particular, (s, 1, 1) configurations can exist in a pure bosonic noncompact case. 

Diamonds can be brought to the bosonic string form with the help of fermionic 
duality transformations. In the 5u(n, m) case this can be done only after introducing 
of a virtual fermionic node. 



3.4 Functional equations for 0l(n|m) spin chain. 



After we gave a description of string configurations that appear in g[(n|m) spin chains, 
it's a time to write down functional equations that involve resolvents for densities of 
these configurations. We will consider first compact su(n|m) case and distinguished 
Kac-Dynkin diagram. The generalizations then can be easily done. 

The distinguished Kac-Dynkin diagram allows, apart usual bosonic strings, only 
for triangle configurations. Let us understand how triangles enter to the string Bethe 
equations for the other nodes. The node (in the notation of (48)) adjacent to 
the fermionic base of s-triangle interacts with triangle as with s-string. The node Vg 
interacts with triangle as with real Bethe root. Bosonic nodes Ufe, 1 < < s, do not 
interact with triangle. Indeed, Bethe equation for these bosonic nodes is written as 



- 1 = 



Qvk-i Qvk+i Qvit 



(57) 



let us represent Q 



Q 



where Qa,s - the Baxter polynomial for the 



s>k 



centers of s-triangles and Q' - the Baxter polynomial for all other Bethe roots. It is 
easy to see that all Qa,s cancels out from (57). 

Procedure of holomorphic projection and introduction of accompanying roots for 
bosonic Bethe roots remains the same as previously. It leads to the following functional 
equations for all purely bosonic string configurations: 



Ca.a' Ra's + Cs.s'Rs 



source term. 



(58) 



This equation has the same form as (47), but identification of the resolvents with string 
configurations is slightly more involved. It should be done in the following way. First, 
assign string configurations to the nodes of the two-dimensional lattice as shown in the 
figure: 



strings of 




(59) 



19 



n s n s 

Figure 3: Figure showing how to derive functional equation for the triangle configurations and 
for the corner. Gray circles show position of real Bethe roots for original Kac-Dynkin diagram. 
Green circles show two auxiliary virtual Bethe roots. Blue region is the place where excitations 
are derived with the help of virtual Bethe roots and therefore not physical. Thus we impose in 
the blue region _R* „ = 0, Rm,s = and i?m+i,n+i = -Rm+i^n+i = 0. We see that after duality 
transformation triangles on the line s = n transform into strings and therefore should obey 
standard functional equations of the form (58). Equation (58) in the corner {m, n} includes 
^m+i.n and and therefore not explicitly expressed through the physical objects. 

String configurations which include fermionic Bethe roots are situated on the line 
s = n. All other nodes correspond to bosonic strings as shown on the figure. If s > n, 
Ra^s is identified with the resolvent for string configurations at {a, s} node and Ra s ~ 
with the resolvent for correspondent holes. For the case s < n holes and particles are 
interchanged: Ra^s is defined to be the resolvent of holes, while R* ^ is defined to be 
the resolvent for particles. 

With the assignments for the resolvents as above, (58) is certainly valid for arbitrary 
a and s ^ n. It is left to show that the same equation holds if s = n, i.e. when it comes 
from the functional equations for triangles. The simplest way to do this is to add the 
virtual fermionic node and to dualize the Kac-Dynkin diagram using inverse of d2 in 
(54). The s-triangles become bosonic s-strings, and therefore the resolvents for them 
should obey the same equation (58). 

The only situation when introduction of virtual fermionic node does not directly 
work, is the longest possible triangle. Indeed, this triangle is situated at position {m, n}. 
It interacts with dualized virtual fermionic node and this interaction is not expressed 
through the physical objects. To express equation in the corner only through the 
physical resolvents, let us consider the following sets of equations (for a,s > 1): 

Caa' Rm-l+a' ,n-l+s + Cgs' Rm-l+a,n-l+s' = ^a,l-Rm-l,n~l+s + Ss,lRm-l+a,n-l- (60) 

oo 

Now, act on these equations with Yl {^s,iD^ + Sa,iD^), the result is 

a,s=l 

oo 

Rm,n + Rm^n = ^ (-Rm-l,n-l+s + Rm-l+s,n~l) ■ (61) 
s=l 

This nonlocal equation is expressed only through the physical resolvents and this is 
the equation that should be considered in the corner. 

Any dualisations of Dynkin diagram will result in reinterpretation of strings as 
triangles and trapezia and vice versa. However, the main equation (58) remains un- 
changed, once this reinterpretation is taken into account. 

We are ready to formulate the main statement of this section. 
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Proposition Consider some Kac-Dynkin diagram for a given g[(n|m) algebra and 

rational integrable spin chain of length L and inhomog entities Xii • • • XL- Suppose that 
all the nodes of spin chain are in the same highest weight irrep with respect to cho- 
sen Kac-Dynkin diagram. If this irrep is such that all Dynkin labels are except one 
which equals to m, then it is possible to formulate string hypothesis. In the thermo- 
dynamic limit Bethe Ansatz equations are approximated by the functional equations of 
the following form: 

^ 1 

Caa'Ra',s + Css'R*a,s' = ^s,l^a,|m| ^ XI ' ^^'^^ 

i=l " 

This equation is most generally valid in the domain called T-hook which is explicitly 
constructed below. In the comers of the T-hook Ra,s o-iT'd Ra,s satisfy equation (61). 

Construction of T-hook and identification of Ra,s and i?* g with resolvents of string 
configurations is done by the following procedure: 

1. First draw Kac-Dynkin diagram on the two-dimensional square lattice by the 
following rule. Put the node with nonzero Dynkin label at position a = 1, s = 0. 
Then for the nodes to the right start moving right and change direction to up or 
right at each node next to the fermionic^^. For the nodes to the left start moving 
left and change direction to up or left at each node next to the fermionic. Note 
that the source node can be bosonic, right or left fermionic or double fermionic 
(noncompact). 

2. Configurations of real Bethe roots are at the position occupied by the nodes of 
the Dynkin diagram. 

3. String configurations with the base node T^^ are on the ray that originates from 
the corresponding Dynkin node T and goes perpendicularly to the link between 
T and adjacent Dynkin node with smaller number of the Bethe roots. 

4. If the lines containing string configurations go horizontally, then Ra^s is the 
resolvent of density of the strings at node {a, s} and i?* ^ is the resolvent for hole 
density. If the lines go vertically, then i?*^^ is the resolvent for string density and 
i?a,s is the resolvent for hole density. 

For main examples of this construction we refer to Fig. 1. Below we add one more 
sophisticated example, where with blue line we marked imbedding of the Kac-Dynkin 
diagram and with red lines we marked collection of string configurations with the same 
base: 

^^Except for the equation (61) for the corner this proposition was initially formulated in [19]. 
Here we add statement about the corner equation and significantly simplify derivation of (62) 
based on duality transformations. 

^^Let us stress that the rule we apply is different from the one in [38], where the path 
of Backlund transformations was encoded by Dynkin diagram which turns at each fermionic 
node. We also apply the rule to turn at the fermionic node, for instance in Fig. 2, but when 
we consider the questions of grading, not organization of string configurations. 

-'^^For simple strings T is just the type of Bethe roots from which strings are made. For 
trapezia, triangles, and diamonds T is the node containing the longest string. 
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(63) 



s 



4 T-hook and all unitary highest weight repre- 
sentations 

In [12] Gromov, Kazakov, and Tsuboi proposed to inscribe a generalization of Young 
tableaux into arbitrary T-hook and gave a rule how to define the unitary highest 
weight representation (UHW) from a given inscription. This rule is distinguished in a 
sense that characters of representations assigned to rectangular Young tableaux satisfy 
bilinear identity 



where characters are nontrivial for {a, s} inside the T-hook domain, equal to 1 on the 
boundaries, and zero outside. 

Quantum generalization of (64) is the Hirota equation on the transfer matrices. 

Here we would like to make a conjecture that the map of [12] from generalized Young 
tableaux to UHW is in fact a bijection. 

The proof of the conjecture is just to take the mapping rule in [12] and to compare 
it with existing classification theorems of UHW in the literature. For su(n), su(n, m), 
su(n|m), and su(n, |m) real forms such proof goes smoothly. Unfortunately, for the 
most interesting su(n, |tn | p) case we found an apparent contradiction in the literature 
(see discussion below) which prevents for the moment to put the conjecture in a form 
of the theorem. Below we will explain what the bijection conjecture should imply and 
for what class of representations the bijection conjecture still has to be clarified. 

Before giving the list of UHW according to the bijection conjecture, we first recall 
the definition of real forms and unitary representations. 

Real form. Consider the Lie algebra q over C. A real form of q is by definition a 
such algebra go over M that 



A possible way to define the real form is to introduce an involutive antilinear automor- 
phism *. It is straightforward to check that generators Jq of the algebra that satisfy 
Jq = —Jq can be used to span go- 
lf instead go is given, one can define * by the property Jq = — Jq and (a J)* = aJ* , 
where o is a complex number, Jo G go, J G g, and a means complex conjugation. 



Xa,s — Xa,s+lXa,s-l + Xa+l,sXa-l,s 



(64) 



= go iXiiR c. 



(65) 
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5u(n) 


5u(n m) 


5u(n, m) 


su(n, m) 


su(n, m p) 


su(n, m p, r) 



Figure 4: List of possible real forms. Cross in a Dynkin node means fermionic 
root, as before. Additional circle around a Dynkin node means noncompact root. 
Number n over the line means n — 1 compact bosonic roots, the same for m, p, t. 

Unitciry representation. The representation T is called unitary for a given real 
form if we can define a hermitian form (positive definite scalar product) on the repre- 
sentation module such that for all J G g 

r(j*) = r(j)+, (66) 

where means hermitian conjugation with respect to the hermitian form. 

Possible real forms. We will focus on the real forms that allow to choose a Borel 
decomposition such that 

E*^ = Eu, E*^ = {-l)<^+<^)E^i, (67) 

where c{i) = or 1; this is a 1^2 grading different in general from the grading The 
root corresponding to Eij is called compact if c(i) + c(j) is even and noncompact if 
c{i) + c{j) is odd. 

Let us consider set of positive simple roots (that correspond to Ei^i^i) together 
with the negative one corresponding to -En+m,i- this set the number of noncompact 
roots is always even, the same is true for the number of fermionic roots. By per- 
forming fermionic duality transformations we can always reduce number of compact 
and fermionic simple roots to at most 2. As a result the full set of real forms and 
corresponding distinguished gradings is summarized in Fig. 4. 

In the following we will consider UHW representations with respect to the standard 
Borel subalgebra defined by one of the cases in Fig. 4^°. 

UHW for su(u). In this case string hypothesis lead to a strip of width n (case (a) 
in Fig. 1). We inscribe a standard Young tableaux in this strip which defines the UHW 

^°There are other choices of Borel subalgebra which lead to a different class of highest 
weight representations. For example, Dynkin diagram of three noncompact bosonic nodes 
still corresponds to the case of su(2, 2) real form. Highest weight irreps in this case are not 
always equivalent to highest or lowest weight irreps defined by Dynkin diagram with standard 
compact-noncompact-compact signature. For some reason such alternative choices of Borel 
subalgebra were not systematically studied in the known to the author literature, though they 
give new instances of UHW [ ] • A class of Bethc Ansatz equations we consider corresponds 
to a standard choice of Borel subalgebra, thus we focus on UHW with respect to this standard 
choice. 
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with weight [mi , . . . , rrin] 





































n 















(68) 



UHW for su(n, m). The maximal compact subalgebra is su(n) ©su(m) ©u(l), with 
the u(l) generator defined by: 

1 " 1 
V n ^ — ' m ^ — ' 

k=l k=l 

We denote the weight by 

[-ai + A, -02 + A, . . . , -On + A, ci, C2, . . . , Cm] (70) 

with condition an > On-i > . . . > ai = and ci > C2 . . . > = 0. The necessary 
condition to get unitarity is that [a„, . . . , oi] and [ci, . . . , Cm] define unitary represen- 
tations for su(n) and su(m) subalgebras respectively. The a— and c— Young tableaux 
are inscribed into a slim hook as follows^^ 



J] 



m 



(71) 



To the a- Young tableaux we added a vertical line of length p which allows map value 
of A as 



^ = -p-Pa-pc, (72) 

where Paipc) is the number of nonzero ak{ck). So, A is minus the total height of the 
figure inscribed into the slim hook. Value of p should be a nonnegative integer if 
P < Po = Max(n — Pa — '\-,va. — Pc — V) and can be arbitrary real if p > pq. 

Comparison with the classification theorem in [10] shows that in a described way 
we parameterize exactly all UHW for the 5u(n, m) real form and standard choice of 
Borel decomposition. 

If we also allow add vertical line to the c-tableaux and define p as a total length of 
a and c vertical lines then restrictions on p are quite natural from the graphical point 
of view. Indeed, if the figure is completely inside the slim hook, it should belong to the 

^^The turning points of the partitions are at the nodes of the lattice, they are shifted from 
the nodes only to improve readability of the picture. 
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integer lattice and the corresponding representation is a short one, analog of the BPS 
state. If the figure can be made, by assigning p only to a or c, such that the distance 
to the boundary is less then one, representation becomes long and p can be deformed 
continuously. 

There is also a nice interpretation of p from the point of view of oscillatory descrip- 
tion. To construct su(n) representation with Young tableaux of height pa, one need 
at least Pa copies of n oscillators to be able performing antisymmetrization. Nonzero 
integer p means that we take more copies of oscillators than needed, namely Pa + p- 



UHW for 5u(n|m). This case corresponds to covariant g[(n|m) representations - 
the ones that can be obtained from tensor power of the fundamental representation. 
All UHW are bijected to the Young tableaux that can be inscribed into a fat hook: 



(73) 



The weight of the representation is given by 

[61, . . . ,6n,Ci, . . . ,Cm], 



(74) 



Note that plays the same role as p + pc in the su(n, m) case. Restrictions on 6„ to 
get UHW are: first, bn > Pc', second, bn should be integer if 6n < vn — 1 or can be 
arbitrary real number otherwise. If b^ > m, one can put Cm to arbitrary nonnegative 
real number. 

UHW^ for su(n-, |m). This case corresponds to contravariant representations of 
g[(n|m) - the ones that can be obtained from tensor power of the antifundamental 
representation. Though we formally have noncompact generators Eij in a sense that 
c{i) + c(j) = 1, all such generators are fermionic and therefore the UHW is finite di- 
mensional. The classification rules are the same as in the previous case, just the role 
of n and m is exchanged and assignment of the weights goes as follows: 



Cl 



m ' 



highest weight is [-Cn, . . . , -ci , ...,-bi] 



(75) 
(76) 
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Figure 5: Possible places of unitarity for 5u(n,|m|p) algebra. 

UHW for 5u(n,|m|p) By fermionic duality transformations the su(n,|m|p) case can 
be mapped to the case su(n, p|m) - the one for the distinguished Dynkin diagram, 
so we won't consider the su(n, p|m) case separately. We choose a grading in which 
the first basis vector Vi is odd: |1| = 1. The parameterization of the highest weight 
representations is the same (up to some relabeling) as in [12]: 
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(77) 



The highest weight is given by 

[A, -02 + A, . . . , -o„ + A, 6i, . . . 6nv, ci, . . . Cp] (78) 

with A = —p — Pa — bi (again minus the total height of the inscription). Furthermore, 
one requirers that p > 0, > pc and that p and b^ arc integers if respectively 
P + Pa^^ — ^oibm<k — 1 (and this corresponds to BPS states), p (b^n) can be 
continuous if p + Pa > n — 1 (6^ > k — 1). Let us translate these requirements to the 
possible values of the fermionic Dynkin labels^^ fi and r2, 



fi f2 




related to the weights by ri = an+ Pa+ P, ^2 = bm + ci. 

If the bosonic Dynkin labels are fixed, the values of fi and r2 that are expected to 
give unitary representations are summarized in Fig. 5. 

The reason that we used "expected to give" is the following. As follows from [41], 
possible places of unitarity shown in Fig. 5 are exactly those that follow from the 

^^here n, m, k mean number plus one of the compact bosonic nodes on the corresponding 
node. 
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first approximation analysis - requirement that action of su(n,|m), su(m|p), su(n, p) 
subalgcbras (separately each) on the highest weight is unitarizable. Subsequent more 
delicate analysis in [41] lead to the conclusion that some of the discrete places of 
Fig. 5 are in fact not unitary. This however contradicts the results of Dobrev and 
Petkova [42] for su(2,|n|2) from which it follows that Fig. 5 describes exactly all UHW. 
For the moment the reason for this contradiction is not clear. 

UHW for 5u(n, m|p, r) According to [41], there is no UHW for the standard choice 
of Borel decomposition. The proof is based on the requirement of simultaneous uni- 
tarity of su(n, m), 5u(p,r), su(m|p), . . . submodules generated from the highest weight 
vector. 

5 Conclusions and outlook 

In this article we derived a set of functional equations (62) which summarizes^^ most 
of our knowledge about thermodynamic limit of unitary integrable rational spin chains 
of gl{n\m) type. We propose to think about these equations as a suitable unification 
point for an analytical analysis of such spin chains, this is also a suitable place to begin 
derivation of the thermodynamic Bethe Ansatz equations, or Y-system supplemented 
with analytical properties of Y-functions^^. Below we list some directions in which our 
work can be useful. 

The set of equations (62), or equivalently (1), is the most symmetric way to formu- 
late the thermodynamic limit of the Bethe Ansatz. For instance, either we study XXX 
spin chain or Gross-Neveu model, we get the same equations (62) with only difference 
in the source term. Therefore (62) make manifest relations between different models. 
And indeed, chiral Gross-Neveu model can be viewed as a special limit of inhomoge- 
neous XXX spin chain. It would be interesting to collect all known continuous limits 
in the literature and classify them using the pattern of (62). 

An unfortunate drawback is that string hypothesis, at least how we formulate it 
here, is not valid for the alternating supersymmetric spin chains [9] . We hope however 
that at least in some reduced form our analysis would be useful also for this case. 

Recent studies of scattering amplitudes [46] showed that it is important to under- 
stand structure of excitations around the GKP string. Identification of these excita- 
tions was done by Basso [47], and as a potential application of (62) wc show in Fig. 6 
how the excitations found in [ '"] are mapped to the nodes of the AdS/CFT T-hook 
pattern. Basically, a particle/hole transformation is made on the node {1,0}. Imme- 
diately from Fig. 6 we see that iS'0(6) isotopic excitations are indeed coming from the 
Bethe Ansatz equations for A-^ Dynkin diagram. Each type of Bethe nodes has one 
source term - the middle node of A3 has {1,0} particles (after particle/hole transfor- 
mation Rifi becomes a resolvent for particles) as a source and left a right nodes have 

^^On the level of equations, not solutions. 
Necessity to solve the AdS/CFT spectral problem at finite volume increased an interest 
in the understanding of the origins of the Y-system. Here we propose an analytic look on the 
problem, in the literature there are also different algebraic approaches [43-45]. 
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Figure 6: Organization of excitations around GKP string described [47] into a 
T-hook domain of [3]. Domain O corresponds to isotopic 5*0(6) excitations, 
domains D^, - to the excitations in the AdS space, fermionic nodes - to 
fermions. Black node is the place were particle/hole transformation is made. 

fermions as a source. For the excitations in the right and left strips of T-hook - it 
is clear that interaction between noncompact excitations goes in the same way as it 
is for XXXii2 spin chain, however interaction with fermions goes through holes, not 
through particles. Therefore, necessity to express holes in terms of particles makes 
Bethe equations for nodes in these strips more complicated. 

Exponentiation property (7) gives us a shortcut way for formulation of the Y- 
system. This observation leads to interesting questions for the AdS/CFT case. The 
AdS/CFT system at finite coupling implies introduction of additional terms in (62) 
for the nodes {1,0}, {1,±1}, and {2, ±2}. Since in the AdS/CFT case the choice of 
monodromy used to define K * piw (1) is not unique, one can get more from the func- 
tional equations (62) than a standard form of the Y-system (8). Presumably, study of 
other monodromies should bring us to results equivalent to [11]. From the other side, 
there is a belief [-3] that the standard form of the Y-system together with simple phys- 
ical requirements is enough to find the finite volume spectrum. If that is true, results 
of [11] should be derivable from physical requirements and therefore finite coupling 
modification of (62) should be derivable as well. This means in particular a new way 
to understand the origin of the dressing factor in the Beisert-Staudacher Bethe Ansatz 
equations. 

A different puzzling feature of the patterns in Fig. 1, derived in this paper from an 
analytical study of the Bethe Ansatz equations, is that these patterns "know" some- 
how about unitary representations of the symmetry algebra. We conjectured that all 
unitary representations can be represented as inscription of certain generalized Young 
diagram inside a T-hook domain. This inscription rule applied for the rectangular 
representations is consistent with bilinear identities on the characters [12] quantum 
version of which are the Hirota equations on the transfer matrices (9). 

It is important to notice that there is no well clarified explanation why the string 
hypothesis gives us knowledge about unitary representations. It is tempting to think 
that string configurations correspond to bound states in a given rectangular represen- 
tation, but explicit realization of this idea is not obtained so far. 

Notice also that we can define physically two different objects: T-functions which 
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are defined through the y-functions by the relation 



a,s+l-'- s,s—l 



(80) 



a,s — 




and the transfer matrices (which are also denoted as T^^^). These physically very 
different objects satisfy however exactly the same equation (9) and the labels {a, s} 

belong to exactly the same domain. In few cases we can make explicit computations to 
show that these objects coincide though without real physical understanding why this 
happens. It would be very important to clarify the reasons for this miracle bijection. 
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A Applicability of string hypothesis 

String hypothesis was always considered as a plausible approximation as it passes the 
check of completeness - explicit counting of states can be shown to give the dimension 
of the Hilbert space of the spin chain model [13,48-50]. However, in literature there 
are known examples when string hypothesis is not satisfied, even qualitatively [51,52]. 
The most severe deviation appears when position u of Bethe roots scales as u ~ L 
with the volume L of the spin chain. This case was in great detail analyzed for XXX 
spin chain in [53]. 

Though the string hypothesis is violated, the numerical studies [54] show that most 
of the Bethe Ansatz solutions obey the string configuration pattern. In this appendix 
we would make a study that supports this point of view. An analytical treatment of 
2-string solution for XXX spin chain is made here, and this gives us a basic idea for 
what regimes string hypothesis can be used. 

We expect that 2-string solution would be a good approximation near origin while 
far from origin deviations from exact string configuration should be severe. The best 
way to track how string approximation depends on the distance to the origin is to 
make moving the center of the string. And this is indeed possible by using the twist 
in the Bethe Ansatz equations. 

Let us add twist (f) to the Bethe Ansatz equations (26) and take logarithm. This 
leads to 



lo "'^ 
2Tri Uk — i/2 2ni 



M 



Uk — Uj+i 
Uk — Uj — i 



I ZTT 



(81) 



where is some integer called mode number. 
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By varying value of (j) we can continuously change location of the Bethe roots and 
in this way change mode numbers Uk of a solution. Large mode numbers correspond 
to the Bethe roots near origin while small mode numbers correspond to large Bethe 
roots. In particular, the point — = corresponds to the Bethe root at infinity. 

We focus now on the case M = 2 and introduce the following notations: 



ui = uq + i/2 + i5, U2 = uq — i/2 — i5. 



(82) 



String configuration can be never exact except at the origin since otherwise log- 
arithmic singularity in (81) appears. Therefore 2-strings separate into two classes: 
2+-strings (with 5 > 0) and 2~-strings (with 6 < 0). Interchange between these two 
classes may happen only at origin. 

Since two Bethe roots are complex conjugated to each other, they should have the 
same mode numbers. Definition of the mode number is unambiguous in the case of 
the 2"'"-strings within standard choice of the logarithm branch cut. For 2~-strings we 
define interaction between Bethe roots using the formal rule log a = Re(log a) if a < 0. 
This makes mode numbers half-integer. 

Consider first the behavior of the 2-string when \S\ is small. Assuming that \6\ is 
much smaller than u we get 



L/2 



4^ 

L 

, 2u 



no <C 1 



, no > 1 



(83) 



So deviation from exact string is exponentially small in L up to no of scale ^fL. 
If no <C 1, equations (81) are approximated by 



27r /L . 
^^0 = I jsign (no) 



(84) 



We see that the maximal possible mode number is ^[-§1, where [x] means integer part 
of X. For even values of L there is an exact string at origin which corresponds to a 
highest vector state of the XXX spin chain. For even values of L 2+- and 2^-strings 
do not interchange while for odd values of L they do. 

Next, consider the behavior of the 2-string 
for large no. When nt — approaches zero Uk 
should become not only large, but also far from 
other Bethe roots. This means that 5 of 2^- 
string should become large. Indeed, assuming 
that no S> 1 and 5 ^> 1 we get 



no 



L-l 
27rnfc — I 



47rni 



24) 



(85) 



2"'"-string significantly deviates from an exact 
2-string at (27rnfc — 0) ~ ^JL — 1, i.e. at scales 
n ~ \/Z or larger. This is consistent with the 
analysis at small values of \5\. 

For 2~-string \5\ cannot become large since we 
cannot pass through 5 = point and negative 




Figure 7: Trajectories traced by 2+ 
(green) and 2^ (red) strings when the 
twist is continuously changed. 
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values of S are bounded by —1/2. The only possibility for the two Bethe roots of 
2~-string to become far from one another is to become two real Bethe roots^^. For 
this they should first collide at the real axis and then separate. After separation we 
no more need prescription for log a with a < 0, this effects into separating of two half- 
integer mode numbers into two subsequent integers. In this way two Bethe roots will 
approach infinity with the interval 2tt of (f). 

In summary, we expect the qualitative dependence of 2-strings on the twist as shown 
in Fig. 7. We support our study with numerical calculations done in Mathematica, the 
code is the following^^: 

BAE[f_]:= Table [((u[k] +1/2) /(u[k] -1/2) )-L Exp [2 Pi I f] == 

-Product [(u[k]-u[i]+I)/(u[k]-u[i] -I), {i,M}], {k,M>] ; 
(* choose one of the initial conditions *) 

(*1*) L=ll;M=2;init={u[l]->0. 457-0. 501 I, u [2] ->0. 457+0. 501 I}; 
(*2*) L=ll;M=2;init={u[l]->0.88 - 0.49 I,u[2]-> 0.88+ 0.49 I}; 
(*3*) L=8; M=3;init=-[u[l]->-0.49- 0.99 l,u[2]-> -0.50, 

u[3]->-0.49+ 0.99 1>; 
step = 1/100; cycle = step Rcinge [1/step] ; ncycles = 2; 
path = (Join 00 Table [cycle+k, {k, 0, ncycles - 1, 1}])+I/100; 
str[path[[l]]] = 

FindRoot [BAE [path [[1]]] .Table [{u[k] , u[k]/.init>, {k, M}]]; 

Do[str[path[[r]]] = FindRoot [BAE [path[[r]]] , Table [ {u[k], 
u[k]/.str[path[[r-l]] ]>,{k, M}] , 
WorkingPrecision -> 15] ; NotebookDelete [temp] ; 
temp = PrintTeinporary[path[[r]] // N] , {r, 2, Length [path] }] ; 

Manipulate [ 

{ListPlot[{Re[#] ,Im[#]}&/@Table[u[k]/.str[path[[r]]] , {k, M}] , 
PlotRange -> {{-4, 4}, {-1.3, 1.3}}, 
PlotStyle -> {Red, PointSize [ . 04] }] , 
Set Accuracy [path [[r]] , 4]}, {r, 1, Length [path] , 1}] 

The main conclusion of the analysis is that string approximation works well for 
Uo < \/L and is inappropriate otherwise. At large values of u mode numbers can be 
approximated by ~ L/(27r«fc). Therefore among all possible complex solutions only 
1 / s/L fraction of them significantly deviates from the exact 2-string configuration. For 
ferromagnetic case however this \j\fL fraction is important since it corresponds to 
the lowest-energy excitations. For antiferromagnetic case the solutions which deviate 
significantly from the exact string are highly-energetic and therefore inessential. We 
conclude that for antiferromagnetic study in the large volume string hypothesis is a 
correct approximation 

We assume that the effects observed on the example of 2-string solution are similar 
in a general situation, in particular approximations (83) and (85) seem not to change 
significantly. This leads to the assumption formulated in section 3.1 as the weak string 
hypothesis. 

^^This is precisely the situation observed in [51]. We see that introduction of twist is an 
effective way to predict this situation. 

^^The twist is made slightly complex to avoid the collision point 
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B A numerical example for string objects 



This appendix aims to give a numerical evidence that more comphcated objects hke 
strings on the nested levels and triangle/trapezium-type configurations do exist. 

First, we consider su(3) spin chain and we will present a string configuration on 
the nested level. The set of the Bethe equations is given by 

BAE := Table [((u[k] +1/2) /(u[k] -1/2) )-L Product [(u [k] -v [i] +1/2) / 
(u[k]-v[i]-I/2) ,-[i,M2}] == -Product [(u[k] -u [i] +1)/ 
(u[k]-u[i]-I), {i,M}], {k, M}]; 
NestedBAE[f_] := Table [Product [ (v [k] -u [i] +1/2) /(v [k] -u [i] -1/2) , 
{i, M}]Exp[2 \[Pi] I f] == -Product [(v[k]-v[i]+ I)/ 
(v[k]-v[i]-I), {i,M2>], {k,M2}]; 
allBAE[f_] := BAE~Join~NestedBAE [f ] ; 

To find an approximative solution we introduce the twist on the nested level and 

choose it to be small and negative. This allows us to choose such Bethe roots on 
the nested level which are large comparatively to the Bethe roots of the momentum- 
carrying node. For the momentum-carrying node we choose a configuration with all 
Bethe roots real. For the nested level we choose a string configuration which can be 
approximated using zeroes of Hermite polynomial [53] : 

L = 108; M = 26; M2 = 5; initf = -1/10; prec = 10; 
Do[n[k]=If [k > M/2, -Round [L/41-M+2 k,Round[L/4]-2 k],{k,M}]; 
BAElog := Table [L Log[(u[k] +I/2)/(u[k] -1/2)] == 

Sum [If [i==k , , Log [ (u [k] -u [i] +1) / (u [k] -u [i] -I) ] ] , {i , 
M}] + 2\[Pi] I n[k] , {k, M}] ; 
initu = FindRoot [BAElog, Table[{u[k], L/(2 \[Pi] n[k])>, {k, M}] , 

WorkingPrecision -> prec] ; 
initv = Table [v[k]->(M+I Sqrt [2M]Root [HermiteH[M2,#]&,k]//N)/ 
(2 \[Pi] (-initf)), {k, M2}] ; 

Then, iterating in twist one can bring the solution to the region where it approaches 
the exact string configuration. 

step = 1/40; path = Table[k, {k, initf, -2, -step}]; 
sol [path [[1]]] = FindRoot [allBAE [initf] , 

Table[{u[k], u[k] /. initu}, {k, M}] ~ Join" 

Table [{v[k], v[k] /. initv}, {k, M2}] , 

WorkingPrecision -> prec] // Chop; 
Do [sol [path [[r]]] = Chop[#, 10" (-prec+4)] &(§ 
FindRoot [allBAE [path [ [r] ] ] , 

Table[{u[k], u[k] /. sol[path[[r- 1]]]}, {k, M}]~Join~ 

Table [{v[k], v[k] /. sol [path [[r- 1]]]}, {k, M2}] , 

WorkingPrecision -> prec] ; 

NotebookDelete [temp] ; temp = PrintTemporary [path [ [r] ] // N] , 
ir, 2, Length [path] }] 
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rep = Rule[a_, b_] :> b; 
Manipulate [{ListPlot [ 

{{Re [#] , Im [#] >&/(§sol [path [ [r] ] ] [ [1 ; ; -M2- 1] ] / . rep , 
{Re[#] ,Im[#]}&/@ sol [path [[r]]] [[-M2; ;-l]]/.rep}, 
PlotRange -> {{-10, 10} , {-2.3, 2.3}>, 
PlotStyle -> {{Red, PointSize[.04]}, {Blue, PointSize[.03]}}] , 
Set Accuracy [path [ [r] ] , 4] } , 
{r, 1, Length [path] , 1}] 

The positions of t;-roots are given by: 

In:= V /O Range [M2] /. sol [path [[-1]]] // N 
Out:= {1.79704 - 1.95446 I, 1.81005- 0.959531 I, 

1.81005+ 0.0404624 I, 1.81004+ 1.04047 I, .79661+ 2.03238 1} 

which is a fairly good string configuration. 

A trapezium or triangle configuration can be simply obtained using a duality trans- 
formation. Let us take the above su(3) example and treat it as a solution of the su(3|l) 
(>*-C)~0~^ ) Bethe Ansatz equations with no Bethe roots at the fermionic node. 
Now, let us perform the duality transformation on the fermionic node which will give 
us Bethe equations for >H3~(2HS) • The Bethe roots of the dualized fermionic node 
are given by: 

In:= Q2[u_] = Product[(u - v[k]), {k, M2}] /. sol [path [ [-1] ]] ; 

u /. NSolve[Q2[u + 1/2] - Q2 [u - 1/2], u] // N 
Out:= {1.7993+ 1.53399 I, 1.79964- 1.45547 I, 1.81004+ 0.540469 I, 
1.81005- 0.459534 1} 

We see that w-roots and the Bethe roots of the dualized fermionic node form a 
(5, 4)-trapezium configuration. The same reasoning can be continued further, and for 
this particular example trapezia configurations up to a (5, l)-trapezium (5-triangle) for 
the su(3|4) algebra can be obtained. 



C Source term in the wing. 

Construction of the main text allows putting the source term in the arbitrary node 
of the vertical strip in T-hook. It would be natural to expect the situation when the 
source term is inside right or left wings. This indeed can be done, however one should 
consider representations with more than one nonzero Dynkin labels. 

To be precise, let us consider su(2, 2|4) case and try to find the situation which 
corresponds to the source term marked by the blue dot in Fig. 8. Expected is to get 
this case from the representation with Dynkin labels (0000032). To be able apply the 
general pattern, we should add two more Dynkin nodes (with zero number of excita- 
tions on them) and perform duality transformation to get the Kac-Dynkin diagram 
with only one nonzero Dynkin labeP^: 




notations here are the same as in Fig. 2. 
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Figure 8: Two different Kac-Dynkin labelings serving the same set of functional 
equations. On tlie riglit tlie situation is tlie same as in main text. On tfie left 
we have two source terms. By demanding absence of excitations at two Dynkin 
nodes depicted by gray on the left, we get su(2, 2|4) T-hook with a source term 
in the right wing. 



It is possible (but not compulsory) to perform duality transformations further 
and get the distinguished Dynkin diagram, the functional equations for which are 
encoded on the right part of Fig. 8. These functional equations obviously contain the 
source term at the desired position. Then, by undoing the duality transformations and 
suppressing auxiliary Dynkin nodes one obtains the su(2,2|4) T-hook. 

The case when the source term is on the boundary of T-hook is also possible. We 
can treat it similarly, by introducing the auxiliary Bethe equations. The only difference 
is that when we rewrite equations in a form that do not contain auxiliary resolvents, 
the source term would not enter (62). Instead, it will appear in the corner equation 
(61). We leave to the reader derivation of how the corner equation should be modified 
in this particular case. 

One may ask, in what situation it is possible to obtain, using duality transforma- 
tions and auxiliary Dynkin nodes, the grading in which only one Dynkin label is not 
zero. To answer this question, let us assume that we've got already this grading and 
that nonzero label is equal to s", positive for simplicity. Then, by applying the rule 
of Fig. 2, we get a lattice with numbers assigned to each edge. These numbers are 
possible weights of En generators'^. We show this lattice in Fig. 9. 

The Dynkin diagrams which can be drawn on this lattice constitute the full set of 
the diagrams reducible to the one with only one nonzero Dynkin label. In particular, 
it is easy to see how to embed the case (86) in Fig. 9 with s' = 1. 



vertical lines correspond to the generators with |i| = 1 and horizontal - with \i\ — 0. 
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Figure 9: Lattice of all possible weights for different choices of the Kac-Dynkin 
diagrams. In the grey regions all the weights are equal to or 's. We draw the 
Kac-Dynkin diagram as a ladder from bottom left to top right turning at the 
fermionic nodes. So, the diagram has nonzero Dynkin labels only when it passes 
through white regions or the center. In particular, maximally only two bosonic 
nodes can have nonzero Dynkin labels, at the places where the diagram crosses 
the borders of the gray regions. 
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